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NOTATION 


a 

E 

g 

H 0 ( ) 

Hit ) 

H 2 ( ) 

h 

k 


m 


nl 


? i 


radius of spherical shell 
modulus of elasticity- 
acceleration of gravity 

/ 3 2 

differential operator defined as Hq ( ) = [ + cot 

\8<f> 2 

differential operator defined as ( ) = (Hq + (1 - 

differential operator defined as H 2 ( ) = {Hq + 2}( ) 

shell thickness 
Eh 3 

12(1 - v 2 ) 

weight of shell per unit surface area 
n(n - 2) . . . 3x1 (n: odd) 

< n(n - 2) . . . 4x2 (n: even) 

, 0 ! ! = (- 1 ) ! ! = 1 


P external force per unit surface area 

P^(cos <j>) , P^ Legendre polynomial of i th order, first kind 

P^ n (cos cf)) , P^ n associated Legendre polynomial of ith order 


u 

w 

a 


in-plane shell displacement 

deflection of shell 

12(1 - v 2 ) — 
h 2 


Y 
e 

V 
v 


viscous damping coefficient 

h 

a 

stress function defined in equations (2) 
Poisson ratio 


)H ) 
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NONLINEAR AXISYMMETRIC FLEXURAL VIBRATION OF SPHERICAL SHELLS* 


Haruo Kunieda 
Ames Research Center 


SUMMARY 


This report presents axisymmetric responses of a nonshallow thin-walled 
spherical shell on the basis of nonlinear bending theory. 

An ordinary differential equation with nonlinearity of quadratic as well 
as cubic terms associated with variable time is derived. The derivation is 
based on the assumption that the deflection mode is the sum of four Legendre 
polynomials, and the Galerkin procedure is applied. The equation is solved by 
asymptotic expansion, and a first approximate solution is adopted. Unstable 
regions of this solution are discussed. 


INTRODUCTION 


The analysis of dynamic behavior is among the most important current 
research bearing on engineering reliability for thin-walled shells. This 
report is concerned with the forced and free axisymmetric nonlinear vibrations 
of nonshallow thin spherical shells. It may be expected that for thin-walled 
spherical shells, especially nonshallow ones, analysis on the basis of nonlin- 
ear theory will reveal such phenomena of response that cannot be provided by 
linear theory. 

The few investigations that have been made concerning this problem of 
spherical shells are limited to the shallow shells (refs. 1 and 2). Reference 
3 has provided the solutions that apply to generalized variational equations 
of motion with the substitution of spatial mode of the deflection derived by 
linear theory. Since the deflection mode of linear vibration appears generally 
as the sum of Legendre functions, with real fractional and complex conjugate 
orders for nonshallow spherical shells, equations cannot be integrated easily 
in variational procedure (ref. 3) . 

In this report the deflection is set as the sum of four Legendre 
polynomials to satisfy the ordinary boundary conditions, and with the Galerkin 
procedure an ordinary nonlinear equation with time as the independent variable 
is derived. 

Although the assumed deflection may not fit the response mode rigorously, 
it appears to be sufficient by comparison with the expansion of the sinusoidal 
function or of the Legendre function with fractional order with Legendre poly- 
nomials. Of course, if the mode of exciting force is the same as that of 
assumed deflection then a rigorous response mode results. 


*This work was carried out while the author was pursuing a National 
Research Council Resident Associateship supported by the National Aeronautics 
and Space Administration. 



The fundamental nonlinear equations and relations for spherical shells 
used here are derived from reference 4 with the assumption of small rotations 
and the neglect of longitudinal inertia terms and nonlinear terms in the 
equilibrium equations of circumferential direction. 

The ordinary differential equation that is finally derived is analogous 
to a nonlinear, single-degree-of-freedom mass spring oscillator. This equa- 
tion, however, has a nonlinearity of quadratic as well as cubic terms of depen 
dent variable, X(t), whereas only a cubic term appears for cylindrical shells 
(refs. 5 and 6). Because of the quadratic term, the solution for X(t) can 
exhibit the nonlinear property of both softening and hardening alternately as 
X(t) changes sign. The range of X(t) over which this behavior occurs depends 
on the deflection mode assumed, on the open angle of the shell, and on the 
thickness-radius ratio. Consequently, it is found that an increase or decreas 
in the amplitude of vibration with an increase in frequency does not necessar- 
ily reflect the effect of a purely softening or hardening spring. 

The solutions of this nonlinear differential equation use the asymptotic 
expansion of this equation, with the expansions in power of thickness-radius 
ratio. The region of the frequency for unstable response is discussed numeri- 
cally as are the effects of thickness-radius ratio, exciting force, boundary 
conditions, and viscous damping. 


GOVERNING EQUATIONS 


The following assumptions may be applicable in the derivation of 
governing equations for thin-walled nonshallow spherical shells. 

1. Large deflection and small strain 

2. Extremely large radius-thickness ratio in comparison with unity 

3. Small rotations 

4. Negligible rotary inertia and transverse shear deformations 
Then, if the angle of rotation is approximated as 

= (u + w°)/a ~ w°/a 

and the longitudinal inertia terms and nonlinear term in equilibrium equations 
in circumferential directions are neglected, the following equation can be 
derived from reference 4 as a governing equation for axisymmetric vibration of 
spherical shells. 

H 2 H 2 (w) - £H 2 eO- yjw 0 °H 2 (T) - |h 2 (w)H 0 (w) + + y OO H 0 W - 2f 00 w 00 

O -I 4 44 

w , r T i r mO m a^ •• ar • a^ ~ 

- — k{H 2 (W) } J + - -p W + Y Y- W - yP = 0 (1) 

The stress resultants N 0 and N are related to w and the Airy-type stress 
function 'V as follows: ^ 
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N„ = 


N, = 


_L [V * 0 + » - |h 2 (w)] 

a 

-y [cot ♦)¥° + ¥ - | H 2 (w)J J 


( 2 ) 


The compatibility condition (in terms of w and is 


H 2 H 1 (*0 - (1 - v) ^H 2 H 2 (w) + aEhH 2 (w) + EhjH 0 (w)w O ° - (w°°) 2 - (w°) 2 ]= 0 


(3) 


Equation (3) is derived from the middle surface strain-displacement relations 

1 f o 1 .0,21 

% ' a L U - " + S 1 J 


e 6 = “ [(cot cj>) (u) - w] 


(4) 


in conjunction with equations (2) 


SOLUTIONS 


In the axisymmetric linear problem, the purely torsional motion of shells 
can be solved (ref. 3) ; then the torsionless mode shape of vibration is consti- 
tuted by three fundamental mode shapes represented as Legendre functions with 
small real fractional order and complex conjugate orders, or two Legendre 
functions with real fractional orders of small and extremely large value and a 
Kegelfunction (ref. 3) . The mode shape represented by Legendre functions with 
real (usually fractional) order determines the mode shape of the overall shell 
vibration as shown in reference 3. In low frequency vibration, the mode shape 
corresponding to small order is dominant; whereas at higher frequencies, the 
mode shape corresponding to large order becomes dominant. These modes can be 
expanded in the series of Legendre polynomials. 

Each shell discussed here has a closed apex and one edge (see sketch (a)). 
For convenience of calculation and to satisfy ordinary boundary conditions, the 
deflection of the following form will be adopted: 


w = W(<f>) • X(t) 


o 

■ ' A ° P U0 * A i P Ml + A 2 p „ 2 + A 3V h • = Zav ■ xtt) 

i=0 1 


(5) 
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where 


P Mi = P y .Ccos <f>) 


Aq a given constant 

(i ^ 0) unknown constant that can be determined with the consideration of 
boundary conditions 

arbitrary integers that specify the particular modes selected 


Stress function 'V can be obtained by substitution of equation (5) into 
equation (3) as follows. 


where 


V = Vi + V 2 


J 

'1 - I £ B„ A t P hxtt) 


1 = 0 


2m 


1 ’l 


'2 = Eh 2 b a P h 2 X 2 (t) (j * 1) 
j =0 J J J 


B. = 

l 


(1 - v){2 - y.^^ + 1)} - a 
(1 - v) - y.(y. + 1) 


( 6 ) 


{2 - j(j + l)Hj(j + 1) - (1 - v)} 


y maximum of y. 

l 

aj see appendix A, equation (A5) 


Longitudinal displacement u can be obtained from the following equation 
deduced from equations ( 2 ) , (4) , and (6) . 


o 


u 



~ — ~ H 2 (w) - Eh aw - 2;- (w°) 2 


(1 + V) vnOO 

Eh a 


- i (wO) 2 


2a 


(1 1 v) 

Eh a 


(»! + ’i'2) 00 
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Then 


u = -Ci 


sj 

eKT^X B i A i p Ji hxtt) 


1 = 0 


2P 

£> + v)b jaj + dj }P j X h 2 X 2 (t) 


( 7 ) 


J=2 


where 


a m 


nr 


x m 


-m 


Lm(m + 1) - 2 2 


] 


d o = U 0 + c 0 ) 


c m see appendix A, equation (A2) 


Equations (5) , (6) , and (7) allow the satisfaction of any ordinary boundary 
conditions. After the boundary conditions are satisfied, the only unknown 
value in the above relations is X(t). Galerkin's procedure applied to equa- 
tion (1) yields an ordinary nonlinear differential equation in X(t). The sub- 
stitution of equations (5) and (6) into equation (1) and integration of 
equation (1) for overall shell surface after multiplication of w sin 4>/ [hX(t)] 
will result in the equation 


X" (x) + 2<X'(x) + X(t) + e6 2 X 2 Ct) + e 2 B 3 X 3 (T) = P (8) 

where the dimensionless independent variable r is defined as t = w Q t, and 
the prime (') means differentiation with x. Other constants are as follows. 


3i = 


X 

-i = n 


^0^ + 1) - 2 

+ 1) - (1 - v) 
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l ^ l 
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i) [y i (p i 


+ 1) 


2 ] 
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32 = 


( '2p 3 

SZ) b j a j A i E 0’ ,, i : »o- 

0 = 0 i=0 


+ e — — (ww° sin <j>| 4* i° cot <p + ~ ( w ) 1 } 

Eh 5 1 L a J, W 


2y 3 


S XI C ) A i 


12(1 - v 2 ) frf . n 

v j =0 1=0 

[Mi/ 2] 


[C^i - 1 ) Bi + 2 - yi(ui + 1 )] E(j ,P i ; 4 > 0 ) 


y] [2(Pi - 2q) + l]B ± ECj ,P± - 2q; 4> 0 ) 
□ = 1 


33 = - ~\ [ww° sin 4> (4*2° cot <j> + 4 , 2)]L =6 

\ 12(1 - v 2 ) Eh 5 9 9 o 


2y 2y 


[j/2] 


*xx c k b.a.|(j - l)E(k,j;4> 0 ) [2(j - 2q) + l]E(k,j - 2q;+ Q D 

j =0 k=o 


3=1 


vJ 

£ 


B 0 = / , A.A.Efy^y.;^) 

i,j=0 


2 _ 


2k = 


6i 


Ehg 


ma 2 


Y 2 ga 2 B ° 
Ehm 3 j 


= -L _ii f ° P AiP si 

®i Eh 2 J "i 


sin <j> d(f> 


Ai=o 


and E(i,j;<i> 0 ) is given in appendix B. 
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Equation (8) contains the nonlinearity of both quadratic and cubic terms, but 
as the coefficients of these nonlinear terms are expected to be very small com- 
pared to the coefficients of linear terms in thin shells, the solutions utilize 
the expansion of this equation in power of the nonlinearity parameter e = k/a. 
Note that 32 > 63, and the sign and magnitude of X(t) determine whether the 
structure acts like a hardening or softening spring. 

Applied surface load is assumed to be harmonic in time and is fixed in a 
mode spatially as in the form 

P = p cos cot = p cos ftx ( 9 ) 

ft = w/oJq 

Then a first approximate solution of equation (8) is obtained 


X(t) = e$o + $1 cos £ + cos + e 2 $3 cos 3 £ 


( 10 ) 


where 


5 = ftT + 0 

*0 = -(7 M 2 + f b 2 2 a 3 ) 


*1 = A + | 8 2 A 2 - | e8 2 pA + iff e 2 6 2 2 A 3 - 8 3 A 3 

$ 2 = t 8 2 A 2 + | 8 2 pA + J B 2 2 A 3 

*3 = + 7 " B3 ) a3 



and amplitude parameter A included in those coefficients can be determined 
from 

e2 (? §3 " t g 2 2 ) A3 + + f 6 2 P - i ° 2 ^ A - P = 0 (11) 

The following relations give the angle of phase difference that is due to the 
viscous damping term k. 

7 

p cos 6 = -$ifi 2 + <!>! + £ 2 3 2 ' 3 > i( 2 $o + $2) + e 2 3 3 ‘l’ 1 3 j 
p sin 0 = 


I 
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Note that the solutions given by equation (10) are not symmetrical; that is, 
nonlinear vibration of spherical shells is physically nonsymmetrical with 
regard to the middle surface. Therefore, the amplitude parameter A tends to 
increase or decrease with frequency depending on the value of 8 2 anc ^ £3* it 
does not necessarily mean a hardening or softening of this structure. Neither 
does the increase or decrease of amplitude of X(t) with frequency necessarily 
mean a hardening or softening structure. 

Equations (8) and (10) indicate that nonlinearity of softening and 
hardening appears alternatively according to the change of sign of X(t) (in 
some region of that value) in relation to the values of B 2 an <i $3* when the 
shell is considered a mass. 


STABILITY OF THE SOLUTIONS IN STEADY STATE 


If small disturbances are added to A and 5 in such forms as 

A A Ax 

A = A + ne 


and 


At 


5 = ? + 5e = 5 + 5i 

the solutions of equation (8) can be given as 

X(t) = e;($o + h) + ( $ 1 + $l) cos 5 + e ($2 + 5 2 )cos 25 

i 

+ e 2 ($3 + $3)cos 35 


( 12 ) 


Since the values of n and 5 can be assumed to be very small and A is con- 
siderably smaller, the high order terms of n and 5 can be negligible against 
n and 5 , and the next approximation will be applicable. 

cos 35i - 1 , sin 35i - 35 i (A" - 0, 5i m = 0) 

Then equation (12) can be rewritten 

XCx) = X 0 (t) + ^o'C-OSi + X 2 (t) (12a) 

where Xq(t) satisfies equation (10) and X^ (t) is of the form 

X^x) = C$o + $1 cos ? + cos 25 + £ 2 $3 COS 35 

Substituting equation (12a) into equation ( 8 ), and neglecting the high order 
terms of Xi and 5 i gives the following equation 
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Xi" + 2 kXi” + Xi + 2eB 2 XoXl + 3e 2 B 3 X 0 2 Xi - p(sin (X 0 ” + kXq'Hi' = 0 


From this equation the following two relations can be deduced. 


-f2 2 $l + 2 kA$i + + 2e 2 B2 + Ml + $1^2 ^ + ^ $2^1^ 

+ 3e 2 B 3 (|- $x 2 + j e 2 $2 2 ^5 1 + p sin 4> - 2Xfi$ 1 5 1 = 0 

2X$i^ + 2 k$^S 1 + p(cos + 2 k$iX5i = 0 

where 

h -ar^ - *iCne XT ) 

To determine n and ? from these two equations, the determinant of the 
coefficients of n and c, must be zero, and X is given as the root of this 
determinant in the following form: 


XX 2 + YX + Z = 0 


where 

x = 4(k 2 + ft 2 ) 

Y = 4k|\ + ft 2 + B 2 e;2 ( 2<i '0 + $2 + + \ <^2) 

+ j B 3 e 2 ^! 2 + j e$ 2 2 )] 

Z = £l - £1 2 + 2 e 2 B 2 ^ 0 + + ~2 ^2 + 2"^ > 2^^ 

+ 3e 2 B 3 ^j $i 2 + j e 2 $2 2 ^J |^1 - ^ 2 + e 2 B2^2$ 0 + $ 2 ) 
+ ^ £ 2 B 3 $i 3 J + 4k 2 ^ 2 

6 - $1 
*1 
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If a root A of this equation has a real part of positive value, the oscilla- 
tion is divergent. This condition occurs if 

Y ^ 0 and Z < 0 ) 

or | (13) 

Y < 0 ) 

Then it is concluded that the oscillation concerning an amplitude parameter A 
that satisfies conditions (13) with associated ft is unstable. 

NUMERICAL RESULTS AND CONCLUSION 

A few numerical calculations will be shown for hinged (w = M<f> = u = 0) 
and clamped boundary conditions. 

Since a part of plane displacement u contains quadratic terms of A^, 
coefficients A^ are determined as solutions of three simultaneous nonlinear 
equations associated with boundary conditions. These quadratic terms with the 
multiplying factor e can, however, be assumed to be smaller than linear ones. 
Therefore, only the linear part of u will be considered here, and A^ can 
be given approximately from linear equations. Numerical results have shown the 
propriety of this assumption. 

Shell With Hinged Edge 

Deflection mode shape depends not only on the selection of orders of 
Legendre polynomials but also on the open angle of shells and boundary condi- 
tions. Deflection shapes will be shown for various sets of y^ with the 
parameter of open angle <Pq in figure 1. 

Case (a) (0,1, 2, 3) 

Case (b) (0,1, 2, 4) for (y 0 , \i lt y 2 , P 3 ) 

Case (c) (0,2, 3, 4) 

The case of $0 = 120° will be discussed here. Figure 2 with regard to W(<j>) 
and U(<J>) shows that the inertia term of u may not be negligible compared 
with that of w in low order mode shape, despite its neglect in this analysis. 
Figure 3 gives relations between amplitude parameter | A | and modified frequency 
n. 1 The case p = 0 means free vibration. Tendency of increase or decrease of 
| A 1 with increasing ft depends only on the sign of the term 
[( 3 / 4 ) 33 ] ” [(5/ 6)8 2 2 ] ( see eq. (11)) and this does not indicate the stiffness 
of hardening or softening of this system as easily recognized from equation (8) . 
Relations between nonlinearity and the inverse of multiplying parameter l/e = a/h 

^ote that ft = w/u > 0 depends on the ratio 81/80 • Values for 8 . (and 
A^) are given in table 1 . 
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will be given in figure 4 in free vibration. The linear state is e = 0, and 
nonlinearity increases rapidly with increasing e. Figure 5 also shows 
| A | - ft relations of p = 0.1 with the parameter of viscous damping k. 

Large values of | A | decrease with increasing k. When three kinds of |a| exist, 
the smallest | A | increases with increasing k. Simultaneously, the stability 
zone will be indicated on figure 5 for the case of k = 0. Oscillation for 
positive A is unstable in the domain bounded by +UL and +LL. Similarly, 
oscillation for negative A is unstable in the domain bounded by -UL and -LL. 
Figure 5(a) reveals that in the range of ft shown in this figure the follow- 
ing cases appear in turn with increasing ft. 


S = stable oscillation 

(2S,1U) (1S,2U) -► (2S,1U) ■+ (IS) 

U = unstable oscillation 


Variation of the instability zone with k will be given in figure 6 . The 
existence of k makes the unstable region of | A | small. Figure 7 shows the 
angle of phase difference according to k in relation with ft. Finally, time- 
dependent function X(t) will be shown in figure 8 for some fixed ft. Note 
that the maximum | A | of an ft does not necessarily correspond to the one of 
three X(t) that has maximum amplitude, and the term cos C is dominant in 
X(t). Also figure 8 indicates that shell surface does not vibrate symmetri- 
cally with respect to midsurface. There is a slight difference in | A | of 
p = 0.1 and p = -0.1 in figure 3, but figure 8 shows that in the final results 
for X(t) only a difference of phase angle appears. 


All 83 are positive, in these examples, 
82 < 0 in (c) . In all cases X(t) 


82 > 0 in cases 
satisfies the condition 

82 


(a) and (b) , and 


X(t) I < 




Then it can be concluded from equation ( 8 ) that the effect of softening 
and hardening appears alternately according to the change of sign of X(t). 




Cases (a) and (b) 

Case ( a) 

X(t) 

positive 

hardening 

softening 

X(t) 

negative 

softening 

hardening 


The deflection shape W( 4 >) in case (a) is similar to that in case (b) but 
has a different magnitude (figs. 2(a) and 2(b)). However, the final results 
for w( 4>,t) differ only slightly. 

Similar figures are shown in figures 9 through 14 for the clamped edge 
condition. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Sept. 24, 1971 
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APPENDIX A 


DETERMINATION OF a-i AND Ci 


The expressions for the stress function ip and the longitudinal 
displacement u contain the terms a^ and (see eqs. (6) and (7)). These 

terms are evaluated as follows. From equation (5), 


sJ 

= X A ± Py .hX(t) 


3 3 


Cw °) 2 = X X A.A.P* P' h 2 X 2 (t) 
i=0 j=0 J M i 


If equation (A2) is set as 


(w°) 2 = ^2 C k P k h 2 X 2 (t) 


the coefficients of this finite series are determined with Legendre polynomials, 
From the expansion (ref. 7) , 


C P 

n,m r r 


it follows by differentiation with respect to <f> that 


P 1 P 1 


■ £ 


D P 

n,m r r 


r= n-m 


i,j D k = 1 {i(i + D + Kj + i) - k(k + l)} i}j c k 
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r 


. .c. 

1,J k 


R i-j+k R j-i+k R i+j-k 
2k + 1 2 2 2 


i + j + k + 1 


R « . * 
1+3 +k 

2 


R = C2s -1)11 
s s ! 


for integer values of s. 


(- 1 )!! = 0 !! = 1 


and if |i - j| is even, then . .D, = 0 (k is odd) 

1,3k 


|i - j is odd, then 


. .D, = 0 (k is even) 
1 > 3 k 


Then the coefficients Cr of equation (A3) can be determined as 

3 3 


= T y A. A. D 

i=0 j=0 J 


At cos (J) = 1, 


(A4) 


P ' = 0, P m = 1 (m is arbitrary integer) 


m 


These coefficients then satisfy the relation 

2 y 2y 2y 

E*- • E 

k-0,2,4,. • • k=l,3,5,. • 


C. =0 
k 


k=0 


(A5) 


The nonlinear terms of compatibility of equation (3) can be expanded in a 
series of Legendre polynomials. 


2y 


H 0 Cw)w°° - (w°°) 2 - (w°) 2 =^a.P. 


(A6) 


1=0 


The left term of equation (A6) can be represented by substitution of equation 
(A3) 
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Then 



2y 

cos <J> p i 

2 sin 4 / , r r 
r=o 


because 


a. 

1 


-C. 

l 


2 i 


-is 


CP 1 ?, cos <|) d<f> 
r r x 



+ 


2 i + 1 
4 


2y 

s 

-r=-i 




i+r+i 


1 } 



cos <|> dc() 


0 


J _ . 2 .i_ 

1 2 i + 1 


1 


(-D 


i+r+l 


r + 1 - i 
r = i 

r > i 


(A7) 


From equations (A 5 ) and (A 7 ) 


a.i - 0 
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APPENDIX B 


DETERMINATION OF EfiJ^o) 


The coefficients appearing in equation (8) depend on the quantity 
E(i,j;<l>o) that is evaluated below. 

/•4>0 

E(i,j;<f>o) = I P^j sin 4> d(f) 

*'0 

i r i^o 

= ^-r. ... 1 srr- sin <KP. X P. - P. 1 ?.) 

1(1 + 1) - 3T3 + 1) L i 3 J 1 Jo 


= ... 1 [sin *(P,- 1 P* - P - X P- , (i / j) 

1(1 + 1) - j (j + 1) L 13 J 1 J 4>=4>o 


A 

i Ci ,i ; 4»0 J = I p i p i sin <t> d<t> 


i + Cl/2) 


i + (1/2) 


[i/2] 


J pd 3 0 <t>0 

P. 4r P . , d4> - ly [p.P. cos 4>1 
I 1 d<j) 1-1 2 L 1 1 Jo 

A 


J - y {P^i COS <j>} 


y; [ 2 (i - 2 P ) + i]E(i,i 


<f>=<)>0 


- 2p;<t> 0 ) 


by using the following expansion 

[i/2] 


aj p i-! - T] (2 W - 2p> * i)Pi. 2p si » 

p^i 


* [i/ 23 : gaussian symbol 
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TABLE 1.- NUMERICAL VALUES FOR CONSTANTS USED IN 


EXAMPLE CALCULATIONS 


(P0»1 J 1»P2»1 J 3) 

Hinged edge 

Clamped 

edge 

(0,1, 2, 3) 

Case (a) 
(0,1, 2, 3) 

Case (b) 
(0, 1,2,4) 

Case (c) 
(0,2, 3, 4) 

A o 

-0.2 

-0.2 

-0.2 

-0.2 

Ai 

-.42379 

-.35341 

.73030 

.035121 

a 2 

-1.19505 

-.87531 

1.85114 

.15909 

a 3 

-.36864 

.29794 

1.79404 

.54274 

&0 

.74510 

.43495 

1.93892 

.00310 

61 

.72470 

.78662 

.94502 

36.67234 

62 

4.35677 

2.78319 

-6.08482 

2.40445 

&3 

4.23661 

2.25442 

54.60953 

1.35567 
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(a) y. = (0,1, 2, 3) 


Figure 2.- Deflection and in-plane displacement shapes for a shell of open 
angle <j>o = 120°; a/h = 100, hinged edge. 






(a) y. = (0,1, 2, 3) 

Figure 3.- Relation between amplitude parameter | A j and modified angular frequency ft with various 
magnitudes of external forces p; <j>o = 120°, a/h = 100, k = 0.05, hinged edge. 
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Figure 4.- Variation of | A | -ft relation with respect to radius-thickness 
ratio in undamped free vibration; p = 0, k = 0, <j> 0 = 120°, hinged edge, 
U, = (0,1, 2, 3). 




(a) y. = (0,1, 2,3) 

Figure 5.- Variation of | A | - fi relation with respect to viscous damping k and instability zone in 
nondamping state (k = 0) . Domains limited by +UL and +LL, and -UL and -LL are unstable regions; 
p = 0.1, 4> 0 = 120°, a/h = 100, hinged edge. 



Figure 5.- Continued. 







Cc) y i = (0,2,3, 4) 
Figure 5.- Concluded 





U) y. = (0,1, 2, 3) 

Figure 6.- Variation of instability zone in relation with viscous damping; 
p = 0.05, <j>o = 120°, a/h = 100, hinged edge. 
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(a) u. = (0,1, 2, 3) 

Figure 7.- Relation between phase-difference and modified angular frequency; p 

a/h = 100, hinged edge. 
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(al) y^ = (0,1, 2, 3); fi = 0.945 (three real roots to eq. (11)). 
Figure 8.- Time -dependent variable; a/h = 100, <j> 0 = 120°, hinged edge. 
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(bl) y i = (0,1, 2, 4); Q. = 0.945 (three real roots to eq. (11)). 

Figure 8.- Continued. 
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<t > 2 -41.6974 

<J> 3 1753.550 
A 7.61078 
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-5.5417 

-29.6193 

-490.7163 

-4.9781 


-23.4022 

-2.784997 

-7.4431 

-72.5795 
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(cl) yu = (0,2,3, 4); 0 = 1.020 (three real roots to eq. (11)). 

Figure 8.- Continued. 
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(a2) y^ = (0,1, 2, 3); 0 = 1.020 (one real root to eq. ( 11 )). 
(g Figure 8.- Continued. 
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(c2) - (0,2,3, 4); ft = 0.945 (one real root to eq. (11)) 


Figure 8.- Concluded 
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Figure 11.- Relation between amplitude parameter 1 A ( and modified angular frequency £1 with various 
magnitudes of external forces p; <t>o = 120°, a/h = 100, k = 0.05, clamped edge, = (0,1, 2, 3). 





Figure 13.- Variation of ] A j - ft relation with respect to viscous damping k and instability zone in 
nondamping state (k = 0) ; p = 0.1, tf> 0 = 120°, a/h = 100, clamped edge, p. = (0,1, 2, 3). 
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(a) y^ - (0,1, 2, 3); fi = 0.945 (three real roots to eq. (11)). 
Figure 14.- Time-dependent variable; a/h = 100, <|> 0 = 120°, clamped edge. 
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